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Weighing the black holes in ultraluminous X-ray sources through 
timing 
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ABSTRACT 

We describe a new method to estimate the mass of black holes in Ultraluminous X-ray Sources 
(ULXs). The method is based on the recently discovered "variability plane", populated by 
Galactic stellar-mass black-hole candidates (BHCs) and supermassive active galactic nuclei 
(AGNs), in the parameter space defined by the black-hole mass, accretion rate and charac- 
teristic frequency. We apply this method to the two ULXs from which low-frequency quasi- 
periodic oscillations have been discovered, M82 X-l and NGC 5408 X-l. For both sources 
we obtain a black-hole mass in the range 100 ~ 1300 Mq, thus providing evidence for these 
two sources to host an intermediate-mass black hole. 

Key words: black hole physics - accretion - X-ray: binaries - X-rays: individual: M82 X-l, 
NGC 5408 X-l 



1 INTRODUCTION 

Ultraluminous X-ray Sources (ULXs) are point-like, off-nuclear 
sources which have luminosities greater than ~ 10 39 erg s _1 , 
in excess of that of a ~ lOM© compact object accreting at the 
Eddington limit. The high luminosity, the very soft thermal com- 
ponent often observed in many sou rces (interpreted a s emission 
from a cool accretion disc, see e.g. IfVliller et"al]|2003l and refer- 
ences th erein; for different i nterpretations of the soft X-ray excess 
see e.g. IStobbart et all|2006L and references th erein) and the vari- 
ability on timescales f rom months to hours dSwartz et alj|2004l ; 
iMucciarelli et al.ll2007l) suggest that (some of) these sources may 
be powered by accretion onto an Intermediate Mass Black Hole 
(IMBH) of 100-1000 Mq. Nevertheless, many of the ULXs prop- 
erties can be explained if the y do not emit iso tropically (geomet- 
rically beamed emission, see lKing et alj|200l[). are do minated by 
emission from a relativistic jet (e.g. | Kaaret et al.ll2003l), are accret- 
ing at a supercritical rate (Begelmaril l2002l ; lEbisawa et alj|2003l ; 
Begelman et al. 2006), or a combination of all these. In this case, 
they may harbor stellar mass black holes and may be similar to 
Galactic black-hole binaries. 

A possible approach to study the nature of ULXs is through 
fast time variability. The analysis of the aperiodic variability in the 
X-ray flux of X-ray binaries is a powerful tool to study the proper- 
ties of the inner regions of the accretion disc around compact ob- 
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jects, since observed frequencies might be linked to specific time 
scales in the accretion disc. The best way to do this is by studying 
the analogies of timing features in ULXs with those in stellar-mass 
black holes. Measurements of the mass of the black hole in ULXs 
can then be derived by applying scaling arguments. 

Low-frequency quasi-periodic oscillations (QPOs) have been 
observed in the X-ray light curves of two ULXs, M82 



X- 1 dStrohmaver & Mushotzkvl 120031) and NGC 5408 X-l 
dStrohmaver et alj " 2007 ). The properties of these QPOs (as, 
e.g., fractional amplitude, variability, coherence, amplitude and 
frequency of the underlying noise, freq uency-flux correlation , 
etc.; for a thorough d iscussion see e.g. Muc ciarelli et al.1 [20061 : 
IStrohmaver et alj|2007l) are reminiscent of those o f the most com- 
mon low- frequency QPO in BHCs (type-C OPO. iRemillard et al.l 
120021; see lCasella et alj|2005l and references therein). This associ- 
ation, assuming an inverse scaling of the QPO frequency with the 
BH mass, allows for an estimate of the black-hole masses in the 
range of a few tens of Mq to 1000 Mq . In order to further constrain 
the mass of the black hole in these two ULXs, several authors ex- 
plore different possible scaling relations. lFeng&Kaaretld2007t) use 
the correlation between power-law luminosity and QPO frequency 
in BHCs to estimate a black-hole mass of around 10 3 Mq for both 
ULXs. 

From the correlation between the spectral photon index 
and the QPO frequency (known to ex i st in BHCs, see e.g. 
IVignarca et al]|2003l) iFiorito & Titarchukl d2004l) derived for the 
BH in M82 X-l a mass of the order of 10 3 Mq. From the same 
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Figure 1. Variability plane for stellar-mass BHCs and AGNs. The lines in- 
dicate v oc M M ~ 2 . The upper line is a fit to soft-state obje cts, the lower 
is a fit to the hard-state BHCs only (from. K ordin g et a l. 2007). 

correlation, Strohmav er et alJ d2007l) derived for the BH mass in 
NGC 5408 X-l a lower limit of 600 M@. From the correlation be- 
tween disc temp erature and QPO fre quency (known in the BHC 
XTE J1550-564. TSobczak et al]|200fj) the same authors derived a 
mass in the range from ~1800 to ~5000 M©. Other possible QPO 
identifications would obviously yield to different mass estimates 
(although the se identifications are more unlikely, for a thor ough 
discussion see lMucciarelli et a02006l ; IStrohmaver et al.ll2007h . 

All the mass measurements/constraints mentioned above 
make use of correlations known to exist in stellar-mass black holes. 
This approach is powerful and physically intriguing, although it 
requires a non-trivial assumption: that the slope of these correla- 
tions, derived for the stellar-mass BHCs, remains the same over 
wide ranges of frequencies, luminosities and - possibly - masses. 
In this Letter we explore a new method to estimate the mass 
of black holes in ULXs. This metho d is based on the correla- 
tion between characte ristic frequencies {Belloni et alj20 02). on the 
"funda mental plane" dMerloni et al.ll2003l) and on the "variability 
plane" dKording et al.l2007D of accreting black holes. The first cor- 
relation is known to hold for neutron stars and stellar-mass black 
holes, hence an assumption is required to extend it to the ULXs 
(see i]2. It . On the other hand, the fundamental plane and the vari- 
ability planes are populated by accreting black holes over a wide 
range of masses, from the stellar-mass BHCs to the supermassive 
black holes in active galactic nuclei (AGNs). 



2 THE VARIABILITY PLANE 

iMc Hardy et ail d2006h has shown that soft-state BHCs and active 
galactic nuclei (AGNs) populate a plane in the parameter space de- 
fined by the black- hole mass, accretion rate and characteristic fre- 
quency m the PDS. Kordin g etal]d2007t) recently showed evidence 
for this variability plane to extend to the hard-states BHCs, with a 
constant offset for the frequencies (see Fig.[T}: 

log vt = log M - 2 log M - 14.7 - 0.9 6 (1) 

where 6 goes from for the s oft state to 1 for the hard state (from 
eq. 6 of iKQrding etaLI d2007l) ). The existence of such a variabil- 
ity plane across different spectral states, if confirmed, would sug- 
gest that accretion on BHs is scale invariant. If this is the case, we 
can thus use the correlation to measure the mass of the black holes 



Figure 2. The QPO centroid frequency (vqpo, which corresponds to 
the QPO observed in M82 X-l and NGC 5408) versus the lower high- 
frequency Lorentzia n (vg , which is the one used for the varia bility plane in 
iKording et alj2007l) in neutron stars and BHCs (adapted from Bello ni et alJ 
|2002|) . The solid line shows the power-law fitted to the data (Eq.|5J, while 
the two dotted lines show the correspondent uncertainty of ±0.5 dex we 
adopt in our estimate of vg (see Sect. 12.31 . 

thought to be hosted in ULXs. The frequencies and the luminosi- 
ties observed in ULXs are somewhat intermediate between those in 
BHCs and in AGNs. Hence it is reasonable and relatively straight- 
forward to put them in the variability plane. The strength of this 
method is that it uses a relation which appears to hold for many 
orders of magnitude in frequency, luminosity and black-hole mass. 
Either ULXs follow the same relation of all the accreting black 
holes, thus yielding constraints on the masses of the black holes 
they host, or they do not, thus demonstrating that their accretion 
flow is somewhat different than in BHCs and AGNs. 



2.1 Estimate of the characteristic frequency 

In order to use the variability plane to measure the mass of a black 
hole, we must correctly identify the characteristic frequencies in 
the PDS of the X-ray light curve of the black hole itself. The rel- 
atively robust identification of the QPOs in M82 X-l and NGC 
5408 X-l with the "type-C" QPO in BHCs (see Sec. £Q, helps 
us in this. The frequency of the type-C QPO is known to correlate 
with other characteristic frequen c ies measured in the PDS of BHCs 
dWiinands & van der Klisl 1 1 9991: [Psaltis et al.l [l999l : iBelloni et all 
120021 : Klei n- Wolt & van der Klis 2007, and references therein). In 
particu lar, the frequency of the component used in Kordi ng" et al.l 
(2007) (the lower high-frequency Lorentzian vg) is found to be a 
factor of ~10 higher than the centroid frequency of the type-C 
QPO (see Fig. [2] an d Eq. O. In Figure ([2} we report data from 
IBelloni et all d2002h and the best fit power-law to them, after re- 
moval of points which do not belong to the plot (because of a dif- 
fe rent and/or dubious i dentification of the two frequencies, see § 4.3 
of lBellonietal] d2002h ). 

If we assume that this correlation holds also for ULXs, we can 
use it to estimate v t for M82 X- 1 and NGC 5408 X- 1 from the QPO 
frequency vqpo ■ 

v e w 12.37 x {vqpo) 1023 (2) 

It must be noted however, that this relation is obtained by extending 
to the ULXs the relation in Fig.l[2]l, which is known to hold only for 
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stellar-mass compact objects. This extensions appears to us reason- 
able, given the fact that the relation appears rather tight over a very 
broad range of frequencies (which includes the frequencies of the 
QPOs observed in M82 X-l) and compact-object masses, but it is 
not certain and will need to be confirmed by the next generation 
of X-ray satellites, which will be able to detect the high-frequency 
component vt in ULXs. 

2.2 Estimate of the accretion rate 

Another key point is of course how to estimate the accretion rate 
from the X-ray luminosity. To do this, we need to know how effi- 
cient is the accretion in ULXs. Many observational evidences sug- 
gest that M82 X-l and NGC 5408 X-l are in a somewhat intermedi- 
ate state (see e.g. lMucciarelli et"al]|2006l) . The main indication for 
this is the high amplitude of the QPOs observed in these sources, 
since BHCs do no t show strong QPOs in their soft state (see e.g. 
Belloni et "ai]|2005h . This means that, a priori, the accretion can be 
either efficient (soft-state like) or inefficient (hard-state like), or, 
most probably, intermediate between the two. We can calculate the 
accretion rate under the two extreme hypotheses, thus deriving a 
range for the mass of the black hole. 

2.2.1 Efficient case 

In case ULXs are efficiently accreting sources (77 ~ 0.1), we can 
measure the accretion rate directly from the bolometric luminosity: 



Where the lower limit is due to the fact that we don't know the 
conversion factor from Lx to Ltoi- Applying Eq. I0 and ([3} to 
the Eq. (0 for the soft-state case (9 = 0), we can thus derive an 
expression for the mass of the black hole: 

los (^) > ^ og (S)"Lk log ^^ 7 - 9 (4) 

2.2.2 Inefficient case 

If the accretion is inefficient, the accretion rate is not linearly 
related to the X-ray luminosity. However, inefficiently accreting 
black holes are known to lie on the so- called "fundamental plane" 
jMerloni et al]2003l ; lFaicke et alj2004h defined by X-ray luminos- 
ity, radio luminosity and mass. This relation holds for many or- 
der of magnitudes in masses, connecting the Galactic stellar-mass 
BHCs to the extragalactic supermassive black holes in AGNs. Un- 
der the reasonable assumption that inefficiently-accreting BHs in 
ULXs would also lie on the fundamental plane (see the Discus- 
sion), we can estimate the accretion rate from the 2-10 keV X-ray 
luminosity: 

M » 0.12 (Lx) ' 5 M 0Ai (—) gs" 1 (5) 

V Tjaccr J 

where r) accr is the total efficiency of the accretion flow, including 
both radiative emission and kine tic energy of any ejected matter 
(from eq. 8 of lKording et al. ! 2006). 

Using again Eq. QJ, this time for the hard-state case (9 — 1), 
and replacing i>t and M with our Eqs.[2]and[5] we derive the ex- 
pression for the mass of the black hole for the inefficient-accretion 
case: 



log M = — log L x - ^ log vqpo - 11.226. (6) 

Eq. and l|6} give the bounds for the mass values for a black 
hole in a ULX with an X-ray luminosity Lx and a characteristic 
frequency vi in the PDS. 

2.3 Uncertainties 

The uncertainties on the mass values given from Eq. l(4j and l[6]l 
mainly come from the intrinsic scatter of the three correlations we 
used to derive the mass estimates. These scatters are larger than the 
observational uncertainties in the measurements of fluxes and QPO 
frequencies. 

From Fig. [2] we see that the full scatter in vt is always smaller 
than half a decade. We can thus associate a very conservative uncer- 
tainty of 0.5 dex (i.e. 1 dex scatter) to the values of vi obtained from 
the measurement of vqpo (Eq. In other words, given a mea- 
surement of log(fQPo), the derived value of log(ue) will have a 
symmetric uncertainty A„ = 0.5. This uncertainty obviously dom- 
inates the one on the direct measurement of vq po (which is of the 
order of ~ 1%). 

The uncertainty on the accretion rate estimate is dominated 
from the scatter in the fundamental plane which, to be conserva- 
tive, we estimate to be of the order of ~0.5 dex (see fig. 4 of 
iKQrding et ai] |2006). which translates in a symmetric uncertainty 
A F = 0.25 dex. 

Since we use the variability plane to estimate the black-hole 
masses, the intrinsic scatter of the plane itself gives of course 
an important c ontrib ution to the uncertainty on the mass values. 
Kordi ng et al.l d2007t) extensively discuss the uncertainties of the 
variability plane, and conclude estimating an overall error on the 
normalization of 0.22 dex. The dependent variable in the variabil- 
ity plane is uiM 2 , which translates in a scatter (full range) of 0. 1 1 
dex on the mass (i.e., a symmetric uncertainty Ak = 0.055 dex). 
An uncertainty of ~ 10% (Al ~ 0.1) can be associated to the 
measurement of the X-ray luminosity. 

These uncertainties are clearly independent, hence we can 
propagate them quadratically through the Eq. 10 and ((6). We do 
not associate any error to the estimate of the accretion rate for the 
efficient case. As we do not convert the X-ray luminosity in bolo- 
metric luminosity, we are considering a lower limit for M. 

The total uncertainties Am on the values of log(M) are thus: 




for the efficient case (Eq.|4), and: 



Am « ^(M + (M + A « + A f = °- 41 (8) 
for the inefficient case (Eq.[6). 

3 THE CASES OF M82 X-l AND NGC 5408 X-l 

Let us now apply these arguments to the two ULXs for which a 
QPO has been discovered: M82 X-l and NGC 5408 X-l. 

3.1 M82 X-l 

The QPO in M82 X-l was observed at 54 mHz and 112 mHz in 
two XMM-Newton observations in 2001 (Strohm aver & Mushotzkvl 
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Table 1. Values of the characteristic frequencies, luminosities and inferred 
masses for M82 X-l and NGC 5408 X-l. 



Source 


v QPO 
(mHz) 


(ergs/s) 


Black-hole mass (Mq) 
Ineff. accr. Eff. accr. 


M82 X-l (2001) 
M82 X-l (2004) 
NGC 5408 X-l 


54±1 
113±2 
20±0.5 


130±13 
170±17 

30±3 


240 +380 

165 +260 
loJ -100 
2 „-+465 


700 +905 
/uu -395 
550+ 710 

570+ 735 



a X 10 ds , in the 2-10 keV range. 

120031) and 2004 dMucciarelli et all 12006k [Pewangan et all 120061) . 
From these frequencies we derive a vg of 0.54 Hz and 1.12 
Hz, respectively (see Sec. ^2. 11 1. After correcting for crowding, 
iFeng & KaareJ J2007t) estimate a 2-10 keV source luminosity of 
1.3 x 10 40 ergs s _1 and 1.7 x 10 40 ergs s" 1 for the 2001 and 2004 
observation, respectively. By inserting these numbers in Eq. l[4j 
we obtain, for the efficient case, a black-hole mass 700^395 Mq 
(2001) and 550±^g Mq (2004). If we assume an inefficient ac- 
cretion, we instead obtain a black-hole mass value of 243tg° Mq 

(2001) and 165±ioo M Q ( 2004 ) from E Q- ©• 

From these values we see that we can put a lower limit for the 
black hole in M82 X- 1 of 95 Mq (smaller values would not be con- 
sistent with the 2001 observation) and an upper limit of 1260 Mq 
(higher values would not be consistent with the 2004 observation). 

3.2 NGC 5408 X-l 

The QPO in NGC 5408 X-l was observed at 20 mHz 
JStrohmaver et alj|2007h . which yields to a vi = 0.20 Hz. We re- 
analyzed the public XMM-Newton data and obtained a 2-10 keV 
unabsorbed luminosity of 3.0 x 10 39 ergs s — 1 . Using these values, 
we derive a black-hole mass of 295l*gQ Mq for the inefficient- 
accretion case and 570lg2o Mq for the efficient-accretion case. 
This translates in a lower limit of 115M© and an upper limit of 
1300A/q for the mass of the black hole in this source. 



4 DISCUSSION 

In Table Q] we report the inferred values for the mass of the black 
holes in M82 X-l and NGC 5408 X-l. It is evident that this method 
supports the identification of both black holes as "intermediate- 
mass black holes". Since we do not know the efficiency of the ac- 
cretion in the two sources, we cannot further constrain the masses. 
However, a comparison with BHCs can give some hints on this. 

QPOs similar to those discovered in these two ULXs are often 
observed in BHCs during their intermediate state. The frequency 
of these QPOs usually increases, and their amplitude decreases, 
as the s ource becomes brighter and the energy spectrum softens 
(see e .g. lCaseila et al.ll2004l ; lBelloni et alfeOOStlHoman & Bellonil 
2005). This softening is usually due to a steepening of the hard 
power law and to an increase of the soft, thermal disc flux. The in- 
crease of the disc flux might be also the origin of the decrease of the 
QPO amplitude, in case the QPO itself arises from the hard power 
law. Hence, the high amplitude of the QPOs observed in the X-ray 
light curve of M82 X-l and NGC 5408 X-l, as well as the pres- 
ence of a very weak disc in their energy spectra, suggest that these 
two sources are on the hard side of the intermediate state. However, 
in case these two sources host an intermediate-mass black hole, a 



relatively low-temperature disc is expected dMiller et alj|200l and 
references therein). 

The QPO in M82 X-l has already been observed at two dif- 
ferent frequencies. Over this small data set the source has shown 
to follow several correlations known to exist in BHCs. Namely, the 
QPO frequency increases and its amplitude decreases as the source 
becomes brighter. This means that during the 2004 observation (see 
TablefJJ the source was probably in a slightly softer, more efficient 
state than in 2001. If we use, for the 2004 observation, the same 
luminosity-to-accretion rate conversion than in 2001, we are actu- 
ally underestimating the accretion rate. This is consistent with the 
fact that we obtain slightly lower ranges of mass values for the 2004 
observation than for the 2001 one. 

In Section §2.2\ we use the fundamental plane of accreting 
BHs to convert the X-ray luminosity in accretion rate for the in- 
efficient case. This method is based on the assumption that ULXs 
lie on the same radio/X-ray correlation as BHCs and AGN do. To 
dat e, radio counterparts have been found only for a few ULXs (see 
e.g. IKording et al1l2005l. and references therein). A full discussion 
about the radio-to-X-ray luminosity ratio of ULXs has already been 
started by other authors, and is beyond the aim of this Letter. Here 
we only note that NGC 5408 X-l is one of th e few ULXs for which 
a steady radio counterpart has been detected. iKording et alj J2005h 
used the values of X-ray and radio luminosity of this source to 
measure the mass of the hosted black hole, obtaining a value of 
~ 1000 Mq, which is roughly consistent with the values we found. 
It must be noted, however, that the observed steep, inverted radio 
spectrum is consistent with a thin synchrotron emission, and not 
with the steady flat radio spectrum observed in hard-state BHCs 
(for a discussion see lSoria et al.l200^) . The slope of the radio emis- 
sion in NGC 5408 X- 1 suggests that the source is on the soft side 
of the intermediate state, when thin synchrotron radio emission is 
often observed in BHCs. This would imply that the mass of the BH 
in this source is closer to the values obtained in the efficient case. 
A bright, variable radio counterpart has been reported also for M82 
X-l. Its strong variability clearly rules out the possibility to use the 
radio flux to place the source on the fundamental plane, since the 
latter is valid only for steady states. The upper limit for the steady 
emission, however, is still consistent with the source lying on the 
plane and a BH of a few hundreds of solar masses. 

4.1 Non-standard accretion flow 

The method described in this Letter is based on the initial assump- 
tion that the accretion flow in ULXs follows the same general rules 
as the one in stellar-mass and supermassive BHs does. If the accre- 
tion onto the black holes in ULXs does not follow general scale- 
invariant relations, the method described here, as well as many 
other scaling arguments, loses validity. Let us thus discuss possi- 
ble alternative scenarios. 

For example, it has been proposed for the ULXs to be sources 
similar to SS433, but with a different viewing angle. In this sce- 
nario, most o f the X-ray lumin osity emitted from a supercritical ac- 
cretion disc (Beselman 2002) w ould be ge ometrically collimated 
(Begel man et al .1 2006: Poutanen et al. 200^). The main concern, 
when invoking a beamed emission (either geometric or relativistic) 
is whether QPOs can be preserved. Any physical interpretation of 
the ULXs involving beaming will need to demonstrate that rela- 
tively coherent oscillations in the X-ray flux do not get smeared out 
by the beaming itself. 

Without applying any beaming correction, i.e. remaining un- 
der the assumption of roughly isotropic emission, there are only 
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two possible ways for a stellar-mass BH ( 20 Mq) to reach 
X-ray luminosities as high as 10 40 ergss -1 : either the accre- 
tion rate is super-Eddington, or the accretion rate is still below 
the Eddington limit, but the accretion flow is extremely efficient, 
as to produce such high luminosities. For example, it has been 
suggested that ULXs host super-Eddington accreting, stellar-mass 
black holes, somewhat si milar to the microquasars in our Galaxy, as 
GRS 1915+105 (see e.g. lKindl2002T) . However, GRS 1915+105 it- 
self lies on many of the relations used here and in literature to mea- 
sure the black-hole mass in ULXs, as well as other very luminous 
microquasars do. The very highly, sometimes super-Eddington ac- 
creting Galactic microquasars appears to follow the same general 
relations of the sub-Eddington accreting black holes. 

To discuss the case of extremely efficient accretion, let us take 
the most extreme case of an efficiency of 50 % (maximally rotating 
Kerr BH) and of a bolometric luminosity equal to the X-ray lumi- 
nosity. Under these extreme assumptions, from Eq. 10 we would 
obtain M = 300± 4 ?° Mq (for the 2004 observation of M82 X-l). 
This demonstrate that, even in case of an extremely efficient accre- 
tion, the method presented in this Letter yields a BH mass higher 
than 130 Mq. 

More observations are clearly needed in order to test the valid- 
ity of the method described here. Detecting more QPOs from M82 
X-l and NGC 5408 X-l, hopefully at different frequencies and/or 
luminosities, will show whether these two sources lie on the vari- 
ability plane or not. The current data strongly argue for M82 X-l 
and NGC 5408 X-l to host intermediate-mass black holes, both 
with a mass from ~ 100 to ~ 1300 Mq . 
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